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Abstract 

In the quantization of a rotating rigid body, a top, one is concerned 
with the Hamiltonian operator La ~ aoLx + ctfLy + o^Li, where ao < 
Q-i < 02. An explicit formula is known for the eigenvalues of La in the 
case of the spherical top (cti = 02 = cts) and symmetrical top (oi = Q2 7^ 
as) [LL] . However, for the asymmetrical top, no such explicit expression 
exists, and the study of the spectrum is much more complex. In this 
paper, we compute the semiclassical density of states for the eigenvalues 
of the family of operators La ~ aoL^+aiLy + a2Ll for any ao < ai < 02. 

1 Introduction 

Let C M.^^ be the 2-sphere and let — A52 be the constant curvature spheri- 
cal Laplacian on 5*^. It is well known that the spectrum of — A52 consists of 
eigenvalues A given by 

Afe = fc(fc + 1), k = 0,1,2,.... 

Moreover, the eigenspace corresponding to is of dimension 2A; + 1 and a 
basis of eigenfunctions is obtained by taking the standard spherical harmonics 
of degree k, i.e. 

rr(^:</>) = ^r(cos0)e™*, |m|<A:, 

where P™ is the associated Legendre function of the first kind. For a more 
detailed treatment of the spectral theory of A52 , we refer the reader to [Foj . 

From the fact that the eigenvalues Afc = k{k + l) of — A52 are of multiplicity 
2 A; + 1, it is easy to see that spectrum of — A52 has clustering. A nice way to 
illustrate this fact is to observe that for any Schwartz function ip on M, 

1 ^ /^/A^^ , /I 



3 = -k 
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in the semi-classical limit k ^ oo (see [M])- Expressions like those appearing 
on the RHS of (|1.1[) are often referred to as a density of states (DOS) (see 
e.g. jTlj ). Together with the mean level spacings and the pairs correlation, the 
DOS represents a useful quantity to measure the spread of the spectrum. 

In this paper we are interested in computing the DOS for \/—La, associated 
with the quantum asymmetric top with Hamiltonian La, where La is given by 

La ■■= [alLl + alLl + ajLl), 

and where 

Lx = -iiydz - zdy), 
Ly ^ -i{zdx ~ xd^), 
Lz = -i(xdy - ydx). 

Here, we assume that a — (ag, al^a^) £ A^, where 

A^ := |a e : < < < a^j 

is the positive Weyl chamber. It is well known that — A52 and —La are com- 
muting, self-adjoint, elliptic operators on L'^{S^) and therefore possess a Hilbert 
basis of joint eigenfunctions - the aforementioned spherical harmonics Yj^ |BT| . 
Moreover, it is easy to verify that their principal symbols are linearly indepen- 
dent in T*(5'2). For these reasons, we say that A52 and La form a quantum 
integrable system on S*^. 

An explicit formula is known for the eigenvalues of La in the case of the 
spherical top (qi = a2 — 03) and symmetrical top (ai — a2 ^ a^) [LLj . 
Although no such explicit formula exists for the eigenvalues of the asymmetri- 
cal top (ai 7^ 0:2 7^ CKs), the spectrum was recently characterized in terms of 
parameters associated with the Lame equation (cf. proposition 2.2 in [T2]). 

For such a system, it is customary to compute the DOS of their joint spec- 
trum (see e.g. |Chl [Co]). Here, we are simply concerned with the density of 
states measures associated to the operators \/—La- In the following, we denote 
by Ek the eigenspace of A52 consisting of spherical harmonics of degree k, i.e. 
Ek = Span{y^ : m = —k, —k + 1, fc}, and by Pk the projection onto Ek- We 
define the DOS measure associated to the operators by 

dpDsix;k,a) := ^ ^ T ^) ^^'^^ 

Ae(T(v-Pfcio) 

where a{y/—PkLa) denotes the spectrum of \/—PkLa- Clearly, a{^/—PkLa) 
consists of the eigenvalues m < of \J —La associated to the spherical 

harmonics of degree fc. Our purpose here is to compute the density of states for 
the measure dpDs{x, k; a) in the semi-classical regime k — > 00. 
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1.1 Main result 

For any given a e A^, let 5 be the function defined on the rectangle [0, tt] x 
[0,7r/2] by 

g{£_,e;a) = (a? - al) (/3cos^ + - l)sin6l) sin^ + ag. 

2 _ 2 

where = Finally, let , 9; a) = max{0, .g(^, 9; a)}. 

"2 "0 

Theorem 1.1. Let be defined as above. Then, we have that 

w-\im dpDsix,k;a) = — / / F{x; 9,^, a) cos 9 d^d9 

fe^oo TT Jo Jo 

where F is a convex combination of delta functions given by 

Fix; 9, a) = (^^ - ^Vg+it 0; a)) + J'^ (^2: - ^Vg+it 0; a) 

The weak limit is taken with respect to Cc(M^). 

The proof of Theorem 11.11 is given in the third section of the paper. In the 
second section, we show how one can separate the variables for the eigenvalue 
problem ~La'>p — AV' and its connection to the Lame equation. In particular, 
we will show how the spectrum of the operators —La can be explicitly computed 
through the Lame equation. 



2 Separation of variables and the Lame equation 

As we mentioned earlier, — A52 and —La are commuting, self-adjoint, elliptic 
operators on L'^{S^), hence they possess a Hilbert basis of joint eigenfunctions 
that form a class of spherical harmonics. Rather than working with the stan- 
dard spherical harmonics Y^™, we introduce a more suitable class of spherical 
harmonics for our purpose, the so-called Lame harmonics |BT| IWW] . 

In terms of the Euclidean coordinates {x, y, z) G M'^, the Lame harmonics of 
degree k are written as 



■tp{x,y,z) = x 




where 7^ G {0,1} and I7I = 71 -f 72 + 73; the value of I7I is chosen so that 
fc — 7 is even. The values of the parameters 9j are determined by the condition 
Arsi/; = 0. A simple computation shows that the 9j^s must satisfy Niven's 
equation 
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Based on Whittaker- Watson [WW] terminology, we say that tp is of the first, 
second, third or fourth species if I7I = 0, I7I = 1, I7I = 2 or I7I = 3 respectively. 
Note that there is no Lame harmonics of the second and fourth species for k even, 
whereas for k odd, there is none of the first and third species. We will see later 
on that there exists respectively k/2 + 1, 3(A; + l)/2, 3fc/2 and (fc — l)/2 linearly 
independent Lame harmonics of the first, second, third and fourth species. In 
particular, for any positive integer k, there exist 2k + 1 linearly independent 
Lame harmonics, hence they form a basis for the space of spherical harmonics. 

2.1 Sphero-Conal coordinates 

In order to describe the Lame harmonics in greater detail, it is useful to introduce 
a different system of coordinates on S'^, namely the sphero-conal coordinates 
|Sp[ IVo| . We denote these by (ui,U2). They are defined for any given positive 
real constants aQ < of < by the zeros of the rational function 

9 9 7 

Riu) — 5- H 5- H 5- 

u — ckq u — a\ u — a2 

where (a;,y, z) €E K^. From the graph of R{u), it is easy to see that < ui < 




Figure 1: The graph of R{u) for fixed values of x,y,z and q^. The at correspond to 
the vertical asymptotes. The intersections with the li-axis are the two roots of R{u) 
corresponding to the values of m. 

The equation R{u) = is invariant under rescaling (x,y,z) i—f (tx,ty,tz), 
so the coordinates (mi,W2) are indeed coordinates on under the assumption 
a;^ + j/^ + = 1- They take their name from the fact that they can be obtained 
by the intersection of the unit sphere with confocal cones. 



5 



The relations between the sphero-conal and EucHdean coordinates are given 

by 



(wi 


-al){u2 






- ao)("i 




(wi 


- al){u2 


-a?) 


(«i 


- a?)("o 


-«?) 




- al)(u2 






- 





In particular, (wi,M2) form an orthogonal system of coordinates on S^. This 
can easily be seen by considering the vectors = {d^iX, d^iV, dmz) for which 

^ ^ _ x"^ y"^ 

- R{u2) 

= 0. 



2.2 Separation of variables 

The great advantage of sphero-conal coordinates over other coordinate systems 
on S'^ is that they allow us to simultaneously separate variables in both of 
the spectral problems for — A52 and —La (see |Sp| ). For example, in these 
coordinates, the Laplace equation —As^^p — k{k + I)?/; takes the form 



U2 — Ui 



d 



-dip 
dui 



= k{k + l)ip 



(2.2) 



where A{ui) — (ui — aQ){ui — a1){ui — a\). One can then separate the variables 
and write ■0(ui,U2) = V'i(''^i)'02(i*2)- Denoting the separation constant by —A, 
it follows directly from (|2.2p that both ipi and ?/'2 are solutions of the same Lame 
equation 



A{x)iP'l{x) + \a' {x)ij',{x) = hk{k + l)x - X)^,{x) 



(i = l,2). (2.3) 



2 ^ ' 4~ 

From the general theory of Lame equation [WW] , it is well known that the 
solutions of (|2.3p are given by the Lame functions 

1^1 (a;) - Mx) = \x~ al[<^'^\x - a\\'^l^\x - al\'^l^<p(x) (2.4) 

where </> is a polynomial of degree (fc — |7|)/2 with 7 chosen as above. Conse- 
quently, the joint eigenfunctions of — A52 and —L^ are given by 



V'(U1,U2) = ]J 1% 
J = l 



v2 172/2 I 



xlV^^^(u,) (2.5) 
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Note that, up to a constant depending only on the a's and the solutions 0j's of 
the Niven's equations, (|2.5p are the Lame harmonics (|2.ip expressed in sphero- 
conal coordinates. 

Based on these observations, we can now compute the eigenvalues of —La- 
Let E be such an eigenvalue; we will show that E = X, the separation constant 
obtained previously. First, we use the fact that 

-{Ll + Ll + Ll)^ = -Ag-^i^ = k{k + 

to deduce that 

{al ~ a\)L\ij + {a\ - a\)Ll^ = {a\k(k + 1) - E)il). 
In terms of sphero-conal coordinates, we can rewrite last equation as 

d ( /—. — - dij} 



U2 - Ui 



{aik{k + l)- E)^. 



Upon separating the variables, ijj{ui,U2) = ipi{ui)'ip2{u2) , we obtain 

A{u,)4>'^{u,) + ^A'iuMui) = - {i = 1,2). (2.6) 

By comparison of (|2.6p with (|2.3p . we conclude that /i = k{k + 1) and E = X 
as desired. All that remains to prove is that we get all the possible eigenvalues 
of —La in this way. This is a consequence of the following result due to Stieltjes 
and Szego (see [Si], §6.3): 

Theorem 2.1. Let Pq, pi, p2 be any three real positive numbers and let ai, a2, 03 
be any three real distinct numbers. There exist exactly m+l distinct real numbers 
u for which the generalized Lame equation 
2 

A{x)y"{x) + ^ pj W{x - a,)y'{x) = {m{m + 1 + \p\)x - v)y{x) (2.7) 

has a polynomial solution y of degree m. Moreover, the m + l polynomial solu- 
tions obtained in this way are linearly independent. 

Replacing the expression of the Lame function V'i given in ()2.4p into (|2.3p . 
one can easily verify that the polynomial 4> of degree (fc — |7|)/2 satisfies the 
generalized Lame equation 

2 



A{x)cj^"{x) + J2 U + I) - "?)'^'(^) 



7( ik-\i\)ik+\j\ + l)x-X + D{a,j)]cP{x), (2.8) 
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where D{a, 7) = {a^ + )72 + (a^ + a^)7i + (a? + a|)7o + 27071 + 27i72ao + 
270720^1- The values taken by = A — D{a,j) in terms of the different values 
of 7 are given in the table below. 



species 


7o,7i,72 


J/ 


1 


70 = 71 = 72 = 


A 


2 


70 = 1,71 = 72 = 

71 = 1,70 72 

72 = 1,70 = 71 = 


X- af- a'i 
X^al-al 
X — — al 


3 


70 = 0, 71 = 72 = 1 

71 = 0, 70 = 72 = 1 

72 = 0, 70 = 71 = 1 


X — 4aQ — a( — 
A - ofg — 4a^ — al 
A — ag — «! — 4q:2 


4 


70 = 71 = 72 = 1 


A-4(a^ + a^ 



Table 1: The values taken by u 

By Stieltjes' result with pi ^ ji + 1/2, we deduce that there are exactly 
(fc — |7|)/2 + 1 distinct value i/ for which (|2.8p has a polynomial solution of 
degree (fc — |7|)/2. In particular, the number of Lame harmonics of degree fc 
and of specie 1 is fc/2 + 1, of species 2 is 3(fc + l)/2, of specie 3 is 3fc/2 and of 
specie 4 is (fc — l)/2. It follows that for any fc e N, there exist 2fc + 1 linearly 
independent Lame harmonics, so they form a Hilbert basis of L'^{S'^). 

Furthermore, for each fc G N, we also obtain 2fc + 1 values of v (multiplicity 
included) to which correspond by Table 1, 2fc + 1 values of A. In other words, the 
eigenvalues of the linearly independent Lame harmonics of degree fc are exactly 
given by the 2fc + 1 values of A. Therefore, we have shown the first part of the 
following theorem. 

Theorem 2.2. Let a = (ao,^!,^!) G A"^, then the spectrum of the operator 
—La is given by all numbers X appearing on the RHS of the Lame equation 
(|2.3|1 . Moreover, the X's corresponding to the Lame harmonics of degree k lie 
within the interval (ao(^ ~ 3)(fc + 1), a|fc(fc + 4) + 4|a|). 

The second part is an immediate consequence of a result due to Van Vleck 
|Va| where he proves that all numbers v corresponding to the polynomial solu- 
tions of degree m of the generalized Lame equation (|2.7p lie inside the interval 
{alm{m + 1 + \p\),a2m{m + 1 + \p\)). It follows from this and (|2.8p that the 
eigenvalues A lie inside the interval 

min{ag(fc-|7|)(fc+|7| + l)-|-D(a,7)} < A < max{a^(fc-|7|)(fc+|7|+l)+D(a, 7)}. 

7 7 

Since 7^ G {0, 1}, it is then easy to see that 



min{a2(fc - |7|)(fc + I7I + 1) + D{a,j)} > al(k - 3)(fc + 1) 
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and 



max{a^(fc- |7l)(fc+ ItI + 1) + 7)} < alk{k + 4) + i\a\ 

7 



from which the conclusion of the theorem foUows. 



3 Proof of Theorem 1.1 

Based on the different species of the eigenvalues, we partition the spectrum of 
—La into four disjoint subsets ai, trl defined by 

erf := {A : A is an eigenvalue of a Lame harmonics of degree k and of species i}. 



For each fc e N, we denote the eigenvalues of \J—La corresponding to the 2fc+ 1 
Lame harmonics of degree k by 

\Am^ < \A^I^ < • • • < \A^H ■ 

Based on the definition of the ct^, we can decompose dpDs{v^ oi) into four 
disjoints sums, i.e. 



2k 



•■rM 



2fc , . 

»=1 A6 




As we mentioned before, when fc is even, only the Lame harmonics of the 
first and third species will contribute to the sum above, whereas only the second 
and fourth species will contribute when fc is odd. Therefore, we can write 




fc even 



The key observation here is that the eigenvalues can be obtained by simply 
regarding the polynomial solution of the generalized Lame equation (|2.8p . More 
precisely, we introduce the sets Zi, i = 1,2,3,4, defined by 



7k ._ 



There exist A G and 7 e {0, 1}^ such that = A - D{a, 7)}. 



Consequently, the four sums above can now be taken over the sets instead 
of af. That is. 
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Moreover, since lyj is compactly supported, we can approximate uniformly ip 
by smooth functions. Without loss of generality, we may therefore assume that 
(/3 satisfies 

since D{a,^) =0(1). The equation (|3.2p easily implies that 

The asymptotic of the sums in RHS of (|3.3p are obtained through the fol- 
lowing lemma. 

Lemma 3.1. Let vq, denote the m + 1 real numbers for which the Lame 

equation 

1 

A{x)v"{x) + pj W{x - a})y'{x) = (m(m + 1 + \p\)x - v)y{x) 

admits a polynomial solution y of degree m. For any if € Cc(K^), we have that 

^_ ^ ^ = i r ^ Ug+{^,e;a)) cose dOdi + O (-\ 

m + \m J TT Jo Jo V y \mj 

where g+{C, d] a) — max{0, {a\ — Uq) (/3sin0cos^ + (/3^ — 1) sin^ 6*) + Oq}, and 

The proof of Lemma lB.ll is rather long and technical, so we prefer to postpone 
it until the end of the present section. With this lemma in hand, we can now 
complete the proof of Theorem 1.1. As a consequence of Lemma l3. 11 we obtain 
for k even. 
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k/2 



— E 



1 

2A:/2 



1 n ^ 

^^'^ \23k/2 



By Lemma IXTl the first sum in tlie brackets of (|3.4p is equal to 



V2 



and the second sum in brackets of (|3.4p is equal to 



TT /.■it/2 



^0 



^ -g+(^,0;a) cos^d^de + O 7 



(3.4) 



(3.5) 



(3.6) 



Combining equations (|3.5p and (|3.6p . we deduce that 



(3.7) 



1 /■''/^ /I 

-/ / F((p;e,6';a)cos6ld6ld^ + C' - 
Jo Jo 



where the function F is defined by 

1 ^l 



r \2 



F{ip; e; a)^-ip[ -git 0;a)]+-ipl -gi^ 0; a) 



3 /3 



4" V2 



The conclusion of Theorem 1.1 for k even then follows from (13. ip and (|3.7p 
Similarly, for k odd, we have that 



2k 



1 

— y 



j=-k 



1 



2k + 1 
1 

2k + 1 



AGcrJ 

E^ 



E^ 



3fc/2 ^ 



3 ^ 



2 3fc/2 



(fc- 



l)/2 E/(2 2^) 



(3.8) 
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As for the case k even, we apply Lemma 13.11 to conclude that (|3.7p holds 
when is a positive odd integer. 

To complete the proof of Theorem 1.1, it remains to prove Lemma |3. II 

3.1 Proof of Lemma 13.11 

According to Theorem 12. II with oq = — 1, oq = and a2 — (3"^ > 0, there exist 
m + 1 real values i7q, i\n for which the generalized Lame equation 

x{x - f3^){x + l)Y"{x) + [pox{x - (3^) + pi{x + l){x - 

P2x{x + 1)] Y'{x) = (TO(m + 1 + \p\)x - i>)Y{x), (3.9) 

admits a polynomial solution Y of degree m. First, we show that for any 



— tE^H^ =-/ / ^m,9;a))coseded^ + 0[- 



(3.10) 

\ iiir I ti /n Ir. \ III I 

where h{^, e;a) ^ (3 sin 6* cos f + {P^ - 1) sin^ 9. 

The starting point in proving ()3.10|) consists of establishing a three-term 
recurrence relation satisfied by the Lame polynomials Y. In particular, this 
will allow us to obtain the eigenvalues of —La as the those of some tridiagonal 
matrix. 

More precisely, we consider a Lame polynomial of degree m of the form 

m 

i=o 

If we replace the expression for Y{x) into the Lame equation p.9p . we obtain 
the following three-term recurrence relation: 

Aj{p,(3)aj + Bj{p,(3)aj+i + Cjip, (3)aj-i = vaj {j = 0, ...,m) (3.11) 

where a_i = 0, Qm+i = 0, and 

'A,ip, 13) = (/32 - -l + p^)- p2j + f3^poj, 
B,{p,f3) = {j + l){j+pi)/3^ (3.12) 
[Q(p,/3)=M-(j-l)(j-2 + |p|). 

Note that, as a result of the above, Aq — = Cq. These relations are more 
conveniently expressed in matrix form. Indeed, if we introduce the tridiagonal 
matrix A ~ [aij), i, j ~ 0, ...,m, given by 

MMl if^=J■-l 

\iz=3 (3.13) 
if * = i + 1, 



12 



then the three-term recurrence relation (|3.1ip imphes that 



where X — (oq, ai, . 



, CLrn 



AX = -X, 

Throughout the rest of the proof, we denote by 



^ the m + 1 eigenvalues of A. Note that the components of the eigen- 
vectors X are exactly the coefficients of the Lame polynomials Y. 

We will divide the rest of the proof into several lemmas. The first one consists 
of computing the trace of the powers A" for any n E N. 

Lemma 3.2. We have that 



In/21 

E 

3=0 



n 

,3, 3, n - 2j 

m / 

1=1 \ 



- 1)"-2J/32J +0(l) (3.14) 



/or any positive integer n. Here, |?t./2] denotes the greatest integer less or equal 
to n/2. 



Proof of Lemma \3 .2[ We decompose ^ as a sum of three matrices, A = L+D+U, 
where D = idiag(0, ^i, Am) and 





/o 





• • 









/o 


Bo 





• • 





\ 




Ci 





•• 
















Bi 


• • 







1 





C2 


•• 




















M 






















Brn- 























• • 


-1 




U 





•• 


Cm 






^0 








• • 





/ 



The trace of A" is then given by the trace of {L + D + Uy\ When we expand 
last expression, the non-commutativity of the matrices L, D and U implies that 
the trace of A" is the sum of 3" terms of the form 



where Mi = L, D, or U. This is unmanageable in its full generality for arbitrary 
n. However, we are interested primarily in the asymptotic information contained 
in the trace, which allows us to make significant simplifications. 

First, we point out that our need to consider stems from the fact that we 
will use polynomials to approximate the continuous function if in Lemma 13.11 
Thus, we need to extract asymptotic information about Tt{A") for fixed, but 
arbitrary n. In our case, we will ultimately be taking a limit m oo for fixed 
n, and so in this limit, n/m — )■ 0. 
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Secondly, we exploit the fact that the terms M1M2 ■ ■ ■ Mn are products of 
matrices, each being lower diagonal (L), diagonal (D), or upper diagonal (JJ). 
This allows us to make definite statements about the zero structure of the matrix 
products, i.e., the entries that are necessarily zero in the matrix product. For 
example, multiplication on the left or right by a diagonal matrix preserves the 
zero structure: LD and DL are both lower diagonal if L is. The analogous 
statement holds for UD and DU. The effect of multiplying by L or t/ is only 
slightly less simple. In fact, as far as the effect on zero structure is concerned, 
L and U behave like quantum mechanical creation and annihilation operators 
(respectively). In detail, if we denote by Rm (resp. Lm) the operation of right 
(resp. left) multiplication by a matrix Af, then for any matrix B: 

(i) RijB corresponds to shifting all columns of B one place to the right: 
coli-|_i(i?(7_B) — coli(-B), creating a zero column in the first colmim. 

(ii) RlB corresponds to shifting all columns of B one place to the left: 
co\i^i{RijB) = co\i{B), creating a zero column in the last column. 

(iii) LijB corresponds to shifting all rows of B up one place: YOWi^i{RijB) = 
row,;(i?), creating a zero row in the last row. 

(iv) i[/_B corresponds to shifting all rows of B down one place: ravii+i{RuB) = 
TOWi(B), creating a zero row in the first row. 

As a result, the diagonal of a term M1M2 ■ ■ ■ Mn in A" will have zero trace 
unless the number of factors j of L is the same as the number of factors of U. 
The remaining n — 2j factors must all be D. Thus, many of the 3" terms do not 
contribute to Tr(^"). 

The last issue concerns the lack of commutativity in the terms that do con- 
tribute to the trace. Some of these terms are of the form 

{LU)W--'^,j=0,...,ln/2l (3.15) 

Since LU and D are diagonal, the trace is particularly simple to compute in the 
case of the canonical terms ()3.15p : 

Tr(MiM2---M„) = TriiLUyn''-^^) 

Noncanonical terms will differ from canonical terms only at order 0{n/m) = 
0{l/m), and so for asymptotic purposes, we may assume that all terms have the 
canonical form (I3.15|) . To see this, note that multiplication of matrices of the 
form L, D, and U constitutes a shifting of their rows and columns. For terms 
with n factors, the number of shifts is at most n. Being products of matrices that 
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are (lower, upper) diagonal, the noncanonical terms will yield sums of products 
of the form 

where Tp, Ag e {Ai/n, Bi/n, Ci/^i \ l — 0, 1, ... , m} and |p — <z| = 0{n). As an 
example, 



4 



/?^(/3^ — 1)-^ T- w " 
/32(/32-l)4 + of- 





^ 1 












f 1 


o{ 






i m 



Since there are exactly ^ "_2j) matrices M\M-2, ■ ■ ■ Af„ that contain j fac- 
tors of X, j factors of U and (n — 2j) factors of D, we finally deduce that 

I"/2l . s 

Tr(A") = y ( " , )Tr((LC/)^i?"-2j)^0(^) 

I"/2l / „ \ ™ / ,-2 \ J / -2 \ »-J 

- g(,,,„_,,)i:0-;^) y ("-ir-"^' 

(3.16) 

This completes the proof of the Lemma. □ 

The next result deals with the inner sum 2_/i=i ~ i^] ( ^ ) i'^ (|3.16p . 
As the next lemma shows, this sum is asymptotically given by a Beta integral. 

Lemma 3.3. We have that 

^ g (i - It)' = 5^ + - ^ + 1/2) + o (i) (3.17) 

where B(p, g) zs t/ie standard Beta integral defined by 



B{p,q)^2 / cos^P-i 61 sin^?-! 9 d9. 
Jo 



Proof of Lemma [5T51 - This is obvious. The LHS of (|3.17p is a Riemann sum for 
the fmiction (1 — x^)^(x^)"^-' on [0, 1], hence 

1 / ,-2 \ J / -2 \ /.I /I 



TO — ' \ TO' / \ TO' / /n \ . 



15 



The conclusion of the lemma follows by making the substitution x = sin 6 and 
using the trigonometric representation of the Beta integral. □ 

As a consequence of (|3.16p and Lemma [3^ it follows that 
-Tr(A") = ^Y.{-) 

-I In/21 . s 

2 ^ \j,j,n~2jj 

+o(-]. (3.18) 



m 

In order to evaluate the sum inside the integral sign, we use the sine function 
defined by 

1 for x = 0, 
^ for X ^ 0. 



sinc(a:) 



The key point here is to observe that sinc(7ra;) = when x is a non-zero integer, 
and that sinc(O) = 1. Using this function, we can then replace the sum in (j3.18p 
by the more appropriate sum over multi-index 7 = (71, 72, 73) such that I7I — n. 
More precisely, we have 



n/2]| 

^ , . " B(j + l,n-j-Hl/2)(/32-l)»-2,/32, 

(3.19) 



E 

- E ('')(/5'-l)^-^/?'^^+^^B(7)sinc(7r(7i-72)). 

where 6(7) := B + ^ + l,n— ^ — ^ + Based on the representation 
of sinc(x) as the integral 

sinc(^(7i - 72)) = ^J e^«(^i-^^) d^, (3.20) 

the RHS of p.l9|) can be written as 

— / (](I3'^ -iV^P'^'+'^^BMe'^^'^'-'^^'' (3.21) 

Replacing B(7) by the expression 

.7r/2 

B(7) ^2 (cos6')''i+T2+i(sin6l)2"-Ti-T=^ d0, 
Jo 
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we can then use the Muhinomial Theorem to evaluate the sum in (I3.2ip . We 
obtain 



|7|=" 



5Z (")('^^ ~ l)''3/3''i+''^(cos6')''i+''^(sin6i)2"-Ti-T2e*«(''i-''^) 



(3.22) 



= (/3 cos f sin 2e+{f3^ -1) sin^ 6y 



If we denote by h{^, 6) = (/3^ cos ^ sin 26* + {0^ — 1) sin^ 9), then equations 
(pTTS]) through ([3?22l) imply that 

— Tr(A") — / / /i"(C, 6*) cos6' ^6* + C — 

The rest of proof of Lemma 13.11 follows by the standard functional calculus 
on the Banach algebra Mm(M), the set of all matrices of order m with real 
entries. However, we can also complete the proof by simply observing that for 
any polynomial P, 

-VTr(P(A))-— / / P{h{i,e))co^d dB di + 0[-\. (3.23) 

Finally, Weierstrass' Theorem implies that for any compactly continuous 
function Lp and any e > 0, there exists a polynomial P with 



svc£,\ip{x) - P{x)\ < e/3. 

X 

\<l>{h{^, 0)) - P{h{^,e))\ COS 61 d9 d^ < e/3. 
From the Spectral Mapping Theorem and (|3.24p we obtain 

< e/3. 



This implies 
2^ 



^Tvi^iA)) - ^Tt{P{A)) 



(3.24) 



(3.25) 



(3.26) 



We choose m big enough in (|3.23p so that 



-VTr(F(A))-— / / P{h{i,e))co^6 dd di 



< e/3. (3.27) 



As a consequence of (|3.25p - p.27p 



1 '"' 1 r'^ 

lim _VTr(^(A)) = — / / ^{h{^, 9)) cos 9 d9 d^ 



(3.28) 
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This completes the proof of (|3.10p for ao = — 1, ai = and 02 = 0^ . In 
the case of interest to us, namely = af^, ai — a\ and a2 — a"^, we use the 
fact that the Lame equation is invariant under affine transformations to make 
the change of variable x ^ x{a\ — a^) + af^. If we let y be the function defined 
by y{x) := Y{x{a\ — a^) + ), then it is not hard to show that y satisfies the 
standard Lame equation 

2 

A{x)y" {x) + ^ - = - '')y^'') 

where v = v{a{ — a^) + af/i. From the fact that — 1 < ^ < /3^, we easily deduce 
that al<^< a\. 

Furthermore, if we introduce the function (^^^(x) := {x{oi\ — ap) + ) for 
any Lp E Cc(M+), then we obtain that 

i=0 i=Q ^ / \ / 

= -E^"(^)+^(-)- (3-29) 
It then follows by (l3?28l) and ([3?29l) that 

iim-V^(^ =— / / Mi,e-a))cose de (3.30) 

where g{£,,9;a) = h{£^,9;a){al — Ofo) + a^. Since p is supported in R+, last 
equation remains valid if we replace 0; a) by g+(<^, 9; a) — max{0, g{^, 9; a)} 
and ip{x) by ip{y/x). This completes the proof of Lemma 3.1. □ 
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